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Abstract 

The effective Hamiltonian for two dimensional quantum wells with rough 

interfaces is formally derived. Two new terms are generated. The first term 

is identified to the local energy level fluctuations, which was introduced phe- 

nomenologically in the literature for interface roughness scattering but is now 

shown to be valid only for an infinite potential well or Hamiltonians with one 

single length scale. The other term is shown to modulate the wavefunction 

and cause fluctuations in the charge density. This will further reduce the 

electron mobility to the magnitude that is close to the experimental result. 
PACS numbers: 68.35.ct, 68.65. +g, 73.20.Dx, 85.30.Vw 
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I. INTRODUCTION 



The vast interest in the physics of charge transport in two-dimensional quantum wells 
stems from potential applications in new devices and subsequent integration with Si-based 
chip technology. Experimentally, it is known that the charge transport inside a quantum 
well is strongly affected by the quality of the well. In particular, it is believed that the 
interface roughness is inherent to the quantum well systems and plays an important role 
for wells at low temperature with small well widthsB. On the theoretical side, starting 
from the seminar work by Kardar, Parisi, and Zhangi, a large effort has been devoted to 
understand the morphology of thin film growth! during the past decade. Nevertheless, 
these works only characterize long wavelength properties of the surface roughness, there 
is no systematic attempt to investigate how the electronic properties, such as the charge 
transport, are affected by the surface roughness. 

The study of the effects of surface roughness on the electronic transport properties has 
a long history, tracing back to the work by Prange and Nee§ on magnetic surface states 
in metals. Later, a more complete model was reconsidered by Andoi. Quite often, these 
works are summarized phenomenologically by introducing a local energy-level fluctuation 
term in potential: (dE/dL) A(r), where E is the energy eigenvalue of the electron, L is the 
averaged well width and A(r) is the local change of quantum well width. Such phenomenol- 
ogy finds its natural application in interpreting the photoluminescence data of GaAs/AlAs 
quantum we iis0i. In this case, it has been established that for temperature less than 80 
K, the linewidth of photoluminescence is mainly determined by the local energy level fluc- 
tuations. Transport properties of two-dimensional (2D) carrier gases at Si/SiC>2 interfaces 
and in semiconductor quantum wells are also shown to be strongly affected by the interface 
roug hness&!§. In these studies, theoretical mobility is also calculated based on the assump- 
tion that the local energy level fluctuation is the dominant effect. It is known, however, 
that the experimentally observed mobility can not be explained solely by the roughness 
roughness. In some parameter regime, one has to introduce, for example, phenomenologi- 
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cally negative impurity charge to account for the extra reduction of the mobility observed 
in experiments0@. 

In this work, we shall systematically investigate the effects of surface roughness. Our 
starting point is an averaged version of the Hamiltonian specialized to the quantum well 
configuration. In an expansion in A(r)/L, two lowest order terms are considered . The first 
term is identified to the local energy level fluctuations (dE/dL)A(r). This term represents 
the mismatch effect of the energy band. It was introduced phenomenologically in the liter- 
ature for interface roughness scattering, but is now shown to be valid only for an infinite 
potential well or Hamiltonians with one single length scale. The other term is shown to 
modulate the wavefunction and cause fluctuations in the charge density. This will further 
reduce the electron mobility. 

The rest of the paper is organized as follows: In Section 2, we lay down the formulation 
of the effective Hamiltonian for two dimensional quantum wells. Some general features of 
the effective Hamiltonian are discussed. In particular, we expand the Hamiltonian to the 
order the A(r)/L and discuss its effects on the basis of single particle states. In Section 3, 
we consider the fluctuation of the charge density caused by the wavefunction modulation, 
and study its consequence on the electron mobility to the order of A(r)/L. In the final 
section, we examine the validity of our approach and conclude. 



II. THEORETICAL FORMULATION 

Let us consider a generic quantum well specified by two interfaces at z — z + (r) and 
z = Z-(r), where r = (x, y) is a two dimensional vector. The average distance between the 
two surfaces is L (see Fig.l). For simplicity, we shall impose the hard wall condition on 
the interfaces. Our formulation is easily generalized to the case when the potential well is 
finite. To investigate the effects that are due to the interface roughness, it is convenient 
to do a transformation that maps z + (r) to L and z_(r) to 0. This transformation is easily 
implemented by 
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/ _ Lz Lz_(t) ,_ 

Z ~A{v) + L A(r) + L' r r ' UJ 

where A(r) = z + (r) — z-(r) — L. After transformation, the wavefunction can be generically 
expressed by \I/ n = ip n (x,y) sin (jjj-z'^) / \J L+A ^ -which satisfies the normalization condition: 
J L+A ^/(ir|* n | 2 = l. 

For typical quantum wells, the Fermi wavelength is about 400A. If L is less than 340A, 
there will be no band crossing at low temperatures, and we can take the average along z 
direction with respect to a given subband, i.e., average with respect to sin ^z' (n will be 
taken to be one). In other words, the more appropriate Hamiltonian to work with is defined 
by 

. _ f L dz' sin tfz'Ksm ?f z' 

After the averaging, we find that an extra potential 5V is introduced to H n = -^-V 2 + 5V 
where the Laplacian is over the (x, y) directions and 

- = a n ^f + Pn 2 VA . {VB _ By A) 



+7n 



A 2 +(vB-jVA"j 2 + 5 n QvA • V + jV 2 A^j . (3) 
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Here A = ^ B = Lz ~ (r) a = (z' 2 -^-\ 8 = (z'-^-\ ^ = <-^-\ and t) = 

iicic J-i — A(r)+L> J-> — A(r)+L' n ~ \ dz' 2 ' n > < n ~ \ dz' 2 ' n ' ^ n ~ \dz' 2 ' n n ~ 

{z'-J^r) n . It is easy to show that S n = —1/2 is generally true. 

If we expand 5V to the linear order of A/L and keep only up to O(z-), we obtain 

5V = E n - — ^A(r) + 6 n — -VA r • V+— V 2 A r) - f3 n V + / , (4) 

L m \L 2L J mL 

where we have identified — 7„/i 2 /2m as E n . Note that the resulting Hamiltonian is invariant 
under reflection: L — z_ — > z + and L — z + — > z_ . Let us first put the last term in the right 
place. For this purpose, we consider the two point correlation function of A(r) 

(A(r)A(r')) = ((z + (r) - L)(z + (v>) - L)) + {z_{v)z_{v')) 
-((z + (r)-L)z4r'))-((z + (v')-L)z4v)). 



If the z + (r) and £_(r) are uncorrelated, the cross terms in the above equation vanish. In 
this case, if one assumes that the surfaces described by z + (r) and z_(r) are statistically the 
same, one obtains that {(z+(r) — L)(z + (r') — L)) = (z-(r)z-(r')) = (A(r)A(r')}/2. In other 
words, both z + (r) —L and z_(r) are of the order of A(r). Therefore, the last term in Eq. (|J) 
is of higher order and can be neglected. In addition to the special case when z + (r) — L and 
z_(r) are correlated, for instance, if either interface is smooth but tilted, i.e., Vz is finite, 
this V(z + — L) ■ Vz_ term will also need to be considered. 

In general, the second term in Eq.@ has no definite relation with E n . When L is the only 
length scale in the Hamiltonian (e.g., an infinite potential well), E n has to be proportional 



to 1/L 2 and — ^ = Eq.(f|) can then be written as 



L dL 

-,2 r 



ivA(r) • V+^V 2 A(r) 



(5) 



Physically it becomes clear that the second term in Eq.(|5|) describes the local energy-level 
fluctuation, which was introduced phenomenologically in the literature for interface rough- 
ness scatteringi'0 but is now formally derived and shown to be valid only for an infinite 
potential well or any other potential with only one length scale. It is easy to check that, to 
the first order in ordinary perturbation theory, the third term in Eq.(|5]) does not contribute 
to the scattering matrix for single-particle states. However, as we shall derive in below, a 
closer investigation shows that this is not correct. Indeed, when this term is combined with 
the kinetic energy, it becomes 

H ' = -*( v+ Ym\ 2 . (6) 

2m\ 2L J K ' 

(a second order term has been neglected). Obviously, its effect is to modulate the wavefunc- 
tion for every particle by 

*(r)-tt(r)exp(-^). (7) 

Note that this result is independent of the depth of the well. There are two consequences: 
First, it causes fluctuations in the charge density. This is a many-particle effect and will 



suppress the mobility. We will analyze it in the next section. Secondly, even for single- 
particle state, it implies that the transmission probability is not one when each particle 
passes by a step in the interface. When combined with the Landaner formula^, it indicates 
that the mobility will be further reduced. We will analyze this effect in the final section. 

In summary, to the first order in A(r), the surface roughness introduces local energy level 
fluctuations and local wavefunction modulation for single particle states. Since the surface 
roughness modulates the single particle wavefunction in a coherent way, the wavefunction 
modulation does not introduce new scattering. In the single particle level, the effect of sur- 
face roughness is entirely contained in the local energy level fluctuations. We shall see in the 
next section that when we include many-particle interactions, the wavefunction modulation 
becomes as important as the energy level fluctuations. 



III. MANY-PARTICLE EFFECT 

In this section, we discuss the many-particle effect that is due to the wavefunction mod- 
ulation. We shall demonstrate its effect on the calculation of electron mobility. The change 
of the wavefunction induces a local modulation in the density of electrons: 



n r. z 



n r 



sin 2 ^') 



cxp( 



(L + A(r))/2 ~~ _Jr v L 

with the understanding that the normalization is done with respect to z. Here n' (r) is the 
2-D electron density after being perturbed by the local energy level fluctuations, A (r). 
It is easy to show that 

2m dE n r „. ^ [ A ( k " P) • el(p_k>r + h.c] 



n' (r) = n < 



7ih 2 k 2 F dL 



k<kf 



k 2 — p 2 



where n = k 2 F /2K is the equilibrium electron density at two dimensions and h.c. denotes 
a Hermitian conjugate of the previous term. We shall assume that the density of positive 
charge background remains unchanged so that the local charge modulation is entirely due 
to electrons. The change of charge density is 
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Sp(r, z) = en(r, z')9(z — z_)9(L + A — z) — en - 



sm' 



L/2 



9 (z) 9(L-z) 



(8) 



For convenience, we shall assume z_ = and neglect the curvature effect due to the roughness 
(for instance, the special case when both interfaces fluctuate while their spacing remains L) . 
To the first order in A(r), the total electric potential 5<f) satisfies 

d 2 _ , snr'f^.:) { . . ,,kz /iik \ A(r) 



sin 2 (—2) 

(V 2 + ^2 )^(r, 2) = 4vrp ind + 47ren ^ 



-2 



1 + — cot 



m dE n 
TT 2 h 2 n dL 



d 2 kj: 



[A(k- p) ■ e ^P- k > r + h. 



k<k f 



k 2 — p 2 



>9(z)9(L-z), (9) 



where /w is the induced charge density. The associated scattering matrix within a given 
subband is given by 

5M(q) = (k\5V\q-k) n 



d 2 rc i Q .rfodzSV(r,z)sm 2 (^z) 



L/2 

2e f L i r±, \ 2, 7171 

— / dzd(p(q,z)sm (—z 
L Jo L 



(10) 



where 8(j)(q, z) is of the order of A. We shall denote J Q L dzS(p(q, z) sin 2 (^z) by <50(q). 

We now express the induced charge density in terms of 5<f)(q). This can be achieved in 
the conventional linear response theory by 

p ind (q,uj = 0,z) = J dz'ReU(q,uj = 0, z, z') e 2 5(f)(q, z), (11) 

where IT(g, u) — 0, z, z') is the polarization insertionlll. If we focus on the n-th subband, the 
one-loop contribution to Il(g, uj = 0, z, z') is 



Re H(q, to — 0, z, z') 



2 ,n7r 16m 



sm {T z )B m l T z, h 2 L2 



V 



d 2 k 



9(1 - k)- 



1 



(12) 



(2tt) 2 v J qk(cos9 + x) 
where q and k are measured in terms of kp, 9(1 — k) is the step function, x = q/2k, and V 
denotes the Cauchy principle value. It is easy to show that 



r 2ir 

V I d9 



o cos 9 + x 







2tt 



if \x\ < 1 
if Ixl > 1 



(13) 



Since the momentum transfer g is always less than 2kp, x < 1 for the range of k integration. 
We find that 



V 



d 2 k 



(2tt) 2 

As a result, we obtain 



9(1 — k)- 



1 



9/ 2 dk 



1 



1 



qk(cos6 + x) Jo 2n q J( q / 2 k) 2 - 1 4?r 



(14) 



Pmd(q,W = o,z) 



Ame 2 . 2 ,riTr ~ 
■— — jSin (— z)50(q) 



(15) 



Substituting the above into Eq.(f|) and performing Fourier transformation on both r and z, 
we find 



(g 2 + k 2 z ) 6<K<h K) = -^^v(k z )S4>(q)-^^A(q) [u (k z ) + v{k z )\ 



+ 



L 2 h 2 

16me dE n 
nLh 2 dL 



d 2 k- 



A(q) 



k<kf 



k 2 — |k — q| 



>v(k z ), 



(16) 



where w (k g ) = ff / L e ifcz ^ sin i^f-z) dz and w(^) = J L e ik * z sm 2 (^fz)dz. The k integration 
has been done in Eq . (|14|) . It is also easy to show that 



ty(q) 



00 c?A; 2 
-oo 2ir 



8(f>(q, k z )v*(k z 



(17) 



Substituting 50(q, fc 2 ) in the above equation, we obtain 

dk z \v(k z )\ 2 



5<f>(q) 



-oo 2tt g 2 + k 2 



167ren A , s 16me 2 
A(q) + 



L 2 



L 2 ft 2 



8me dE n 

q, -i^^r A(q) 



00 rf/c 2 u*(A; z ) • u(k z ) 167ren 
-oo 27 g 2 + k 2 U~ 

The A; z integration can be analytically solved as: 



A(q) 



(18) 



/(») - / 



1 



1 



dfc 2 Hfc z )| 2 

2tt g 2 + fc 2 ~ 2 [g 2 ' 2[g 2 + (2mr/L) 2 ]_ 
dk z v*u 
2tt g 2 + k 2 



+ 



(2nvr/L) z 



-qL _ j 



4g 3 [g 2 + (2n%/L) 



21 2 



(19) 



/n7r 



n7re 



-qL 



nn 



\LJ [gL [g 2 + (2mr/L) 2 ] 3 g 2 [g 2 + (2mr/L) 2 ] 2 q 2 \2mx 
With a little rearrangement, Eq.flTBD gives 



(20) 
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For a narrow quantum well satisfying qL < 2kpL <C 1 (this requires L <C 33A for n 

g and J(q) « 



2 x 10 15 m 2 in quantum wells), I(q) ~ |^ and J(g) ~ — and Eq.(|21|) reduces to the 



standard 2D screening form!: 



ty(q) = [-A(q)] (22) 



(the second term in the denominator of Eq.(|2T|) does not exist in pure two dimensions) where 
q s = 2me 2 /h 2 ~ 1/ (0.25 A) and a = iren . 

However, if qL « 1, one shall have to use the full expression of I(q) and J(q). Since 
— 10 30 m -3 and I{q) ~ 10 _27 m 3 , the second term in the denominator of Eq.(^Tj) 
dominates and 50(q) ~ 2 y^A(q). The resulting scattering matrix within a given subband 



is thus given by 



war - vmr> - £ ( t - [ T + "S^ * m ' sw 



2 



"n( 2 ^l-^^J (23) 

where K is the normalization and 5(g) is the power spectrum of A(q), given0 by (|A(q)| 2 ). 
Given the scattering matrix, we can calculate the relaxation time via the relation 

7jk) = 2^ / dV l M ( k " k, )| 2 ( X - cos$)5(£(k) - £(k')), (24) 

where <3> denotes the angle between the initial and final wavevectors k and k'. The mobility 
of the electron can then be solved by 

P J 4nk B T cosh 2 ((E -E F )/2k B T) v ; 

We see that the interparticle interaction reduces the electron mobility estimated by the 
energy level fluctuations by at least three quarters. Since < 0, the second term in 
Eq. (P3"|) due to wavefunction modulation increases the scattering matrix and further reduces 
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the mobility. The overall reduction of the mobility in comparison to previous approach is 
about 1/4. 

We note in passing that in general, in addition to the above Coulomb interaction, the 
density modulation induced by the surface roughness also affects any interactions that de- 
pend on the electron density. If in the absence of surface roughness, the interaction is 
described by / dr J <ir'n(r)Vo(r, r')h(r'), then formally the effect of surface roughness can be 
simply included by replacing Vo by 

V(r,0»V„(r,0(l-^±^M). (26) 



IV. DISCUSSIONS AND CONCLUSIONS 

In this section, we examine the validity of our approach by studying a simplified rough 
interface: z + = L ■ 9(—x) + (L + A) • 6{x) and z_ = 0, i.e., a step at x = (see Fig. 2). We 
shall directly investigate the solution without using the effective Hamiltonian Eq.(|3]) and 
show that they are consistent with each other. Let us first fix the boundary condition by 
requiring the wavefunction to be travelling waves when far from the step: i.e., *&(x, z) — > 
(e lkx + TZe~ lkx ) sm(irz/L) as x — > — oo and ty(x,z) — > Te ipx smijiz/ (L + A)) as x — > oo. 
When we approach x — 0, higher subbands begin to participate, and the most general form 
for the wavefunction is 

{oo 
(e ikx + TZe~ ikx ) sin f + £ b n e knX sin for x > 

oo n=2 ' (27) 

Te ipx sin ^ + E a { e~^ x sin for x < 

where k n and pi are positive and satisfy the energy conservation condition: 

-^(t) 2 -^(i^) 2 - 2 + (i) 2 - 2 + (i^) 2 - <*> 

The fact that 7r/L is comparable to the Fermi momentum makes these higher subbands 
(with n, I greater than one) correspond to decaying modes when A C L. Also in this limit, 
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since TZ, b n and a/ vanish when A = 0, we expect them to be no greater than the order of 
O(AfL). 

To find the coefficients in ^f(x, z), we apply the matching conditions at x = 

for < z < L 

~ (29) 




d x V{0 + , z) = d x ty{0~, z)for0<z< L. (30) 
By using the completeness relation of sin(7r,s/(L + A)) and Eq. (p9"D , we obtain 

oo 

a, = (1 + K)Y{1, + b n Y(n, I), (31) 



n=2 



where a± = T and Y(n, I) is defined by 



Yin, t) = / az sm sin 

L + AJo L 

s . (32) 



L + AJo L L + A 

2n(-l)» sin ^ 



This implies that Y(n,l) = 0(A/L) when n ^ I, while Y(n, /) = 0(1) when n — l. This 
fact, when combined with Eq.(^TJ), yields 



T = (l + ft)Y(l,l) + 0(A 2 /L 2 ). (33) 

Similarly, using the completeness relation of sin(7T2;/L) and Eq. (|3"0"|) , we obtain 

r i a / oo \ 

ik(l -11) = — ^- [ipTY(l, 1) - 5>jPjY(1,Z)J , (34) 

which, to first order in A/L, reduces to 

fc(l - ft) = £±^pTY(l, 1) + 0(A 2 /L 2 ). (35) 

Lj 

By combining Eqs.(|33"D and (p5|), the coefficient T can be determined as — ^ (l + ifp-)- 
Note that the second term — ^ reproduces precisely the rescaling of wavefunction in Eq.(^), 
while the third term — 57375 is nothing but the same reduction^ of transmitted amplitude 



2k 2 L 3 

in the 1-D Schrodinger equation with a potential barrier of height, §f A (the energy level 



fluctuations in Eq.(|5|)): 
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Although A is assumed to be positive in the above derivations, we have checked that our 
conclusions remains valid when A is negative. 

We have thus seen that the validity of our effective Hamiltonian approach has been fully 
checked. We now briefly re-examine the effect due to the change of the single-particle state 
in the ballistic regime. For a single step as we considered in the above, according to the 
Landauer formula0, the conductance due to the step is given by 

e 2 T 

G = ^ 

7m R 

where T and R are transmission and reflection probability. A simple analysis shows that 
the mobility is given by 
lei T \e\ 



irTirin R nhric 



l(i + 7i 2 /e F L 2 ) 1 

A 



4836 



L(l + n 2 /k 2 F L 2 ) i 
A 



cm 2 
Vsec 



(38) 



When the system has many steps, we simply replace A by ^i- Since A is at the order of 
3-4 A, we can estimate A; by the same order. Therefore, \x is about 10 4 — 10 5 cm 2 /V sec for 
L = 100 A. This number when combined with the contribution from Eq.(|2"3"|) (~10 5 cm 2 /Vsec) 
predicts that the mobility is at the order of 10 4 cm 2 /Vsec, in close to experimental result^. 

In conclusion, we have derived an effective Hamiltonian for two dimensional quantum 
wells with rough interfaces. It is checked to give consistent results for an exactly solvable 
model. Two new terms are generated. The first term is identified to the local energy level 
fluctuations, which was introduced phenomenologically in the literature but the previous 
form is now shown to be valid only when the Hamiltonian has one single length scale. The 
effect of this term on the electron mobility has been discussed before. The other term is 
a new finding, which is shown to modulate the wavefunction and cause fluctuations in the 
charge density. We discuss its effects on the reduction of the electron mobility both at 
the level of the single-particle state and by including the many-particle interactions. An 
estimate of the electron mobility is made and gives rises to correct order in comparison to 
experimental data. 
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Figure Captions: 

Fig. 1 A schematic plot of a 2D quantum well with rough interfaces 
Fig. 2 A simple step at x = can be solved asymptotically. 
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